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The study presented in this paper shows that accounting for the relativistic structure of 
the deuteron allows to explain the ratio of the Drell-Yan pair production cross-section at 
the low Bjorken x off the deuteron and the proton. Thus, the sea quark distributions in 
the nucleon should be studied with accounting for the effects of the relativistic structure 
of the deuteron. The suggested approach reduces theoretical uncertainty in extracting 
the ratio u/d from the data and it is important for the clarification of the nature of the 
sea quark asymmetry in the nucleon. 
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1. Introduction 

The Drell-Yan pair production remains up to now the theoretically cleanest way to 
access sea-quark distributions in hadrons. Necessary stage of such analysis is the 
measurement both with protons and neutrons, in which sea-quarks with different 
flavors are probed. As the construction of a target containing free neutrons is not 
feasible, nuclear targets have to be employed. In the experiment that used the 
simplest nuclear target — the deuteron, the flavor asymmetry of the sea-quark 
distributions in the nucleon has been observed^. 

At the same time, analysis of the deep inelastic scattering off different nuclei 
shows that nucleon valence quark structure changes beyond mass- shell El. Recent 
studies within the Bethe-Salpeter formalism point to the relativistic nature of this 
phenomenon. Relativistic effects should be manifested as modification of the sea 
quark distribution in the bound nucleon and, therefore, will change the contri- 
bution of the sea quark component in the nucleon structure extracted from the 
deuteron data. This, in turn, affects estimations of the anti-quark component made 
by analyzing the deuteron and proton data. 

To clarify the role of the relativistic effects in the Drell-Yan (DY) pair produc- 
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tion, we study this process within the Bethe-Salpeter approach. 
2. Drell-Yan Crossection 

Since the amplitude for the Drell-Yan scattering is defined by the imaginary part 
of the amplitude for the forward scattering (A\T(J^J V )\A) in the cross-channel, we 
can obtain it within the Bethe-Salpeter formalism developed for DIS Within this 
approach the matrix element is defined by the nucleon space-time distribution and 
vacuum average of the T-product of the nucleon fields and nucleon electromagnetic 
current ( see Ref. ^ . Following this method we get the following expression of the 
cross-section for the Drell-Yan scattering off deuteron: 

{ ' y; J (2 7 r) 4 (p2_ TO 2)2^p_ :p) 2 _ m2) - w 

This expression gives the cross-section of the DY pair production off the deuteron in 
terms of the off-mass-shell DY nucleon cross-section <r pN and the nucleon distribu- 
tion function / N / D (P, p). This distribution function is defined by the Bethe-Salpeter 
amplitudes, and together with the denominator it composes the four-dimensional 
momentum distribution of the struck nucleon inside the deuteron carrying the to- 
tal momentum P — (Md,P). In that way Eq. |T]) expresses the nucleon blurring 
that results from the four-dimensional distribution of the nucleon inside the nu- 
cleus. By analogy with the non-relativistic 3D momentum distribution we call it 
four-dimensional Fermi motion. 

Due to the four-dimensional integration in Eq. (TT]) actual calculations require 
information about nucleon cross-section cr pN in the kinematical region of the off- 
mass-shell values of the nucleon energy po (pg ^ p 2 + m 2 ). The off-mass-shell 
behavior of cr pN is unobservable, since then explicit microscopic calculations of the 
cross-section have no experimental reference and strongly model dependent. Thus, 
Eq. (TIJ has to be rewritten in terms of measurable quantities such as the cross- 
section for DY scattering off the physical nucleon. The simplest solution of this 
problem is provided by the integration in Eq. {T]) with respect to po- Analytical 
properties of the integrand in Eq. (TT]) give a way to do it explicitly. Within the 
assumption that a pN (p, q) and f N (P,p) are regular with respect topo the integrand 
contains a second order pole corresponding to the struck nucleon and a first order 
pole corresponding to the spectator. These poles lie in the different half-planes of 
the complex plane po- So, we can choose one of the singularities to perform the 
integration. To express the nuclear hadron tensor in terms of the nucleon structure 
functions it is necessary to choose the second order pole (po — Sn) that corresponds 
to the struck nucleoli. The result of the contour integration in vicinity of the second 
order pole in Eq. |T|) is defined by the derivative of the pole residue with respect to 
Po at the point po = Ej$. 

Doing the integration and using the relation in the Bjorken limit cr(P p , Pd, q) = 
a(xi,X2) we get the following expression for the cross section of the DY pair pro- 
duction off the deuteron: 
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a^(x 1 ,xf) = 



d 3 p 
(2^ 



a ( Xl ,x 2 ) E ^x 2 ^ 



f N / D (M D ,p) 
8M D E^ 2 ' 



where x 2 = xJ?m/(i?N —Pz) is the Bjorken x of the struck nucleon, — m 2 + p 2 
is the nucleon on-shell energy, A^ = Md — 2E^ is equivalent to the total energy 
shift of the struck nucleon due to the binding and Fermi motion. 

The function / n ' d (Md,p) together with the denominator composes the three 
dimensional momentum distribution of the struck nucleon inside the deuteron. Ac- 
cording to the normalization condition for the Bethe-Salpeter vertex function this 
distribution satisfies the baryon and momentum sum rules - and coincides with 
the usual nuclear momentum distribution. Thus, the first term in Eq. I[2|), which 
results from the derivative of the propagator of the nuclear spectator, expresses 
contribution of the conventional 3D-Fermi motion. 

The second term with da pN (o^nV'^n results from the nucleon DIS amplitude 
derivative: 

da$(p,q) _ dx* da^jxux™) 
dpo dpo dx 2 

This term is responsible for the deviation from unity of the ratio of the cross sections. 
Its contribution is proportional to the coefficient A^; which, therefore, characterizes 
contribution of the four-dimensional Fermi-motion. 

It is important to note that the struck nucleon pole gives the factor 1/A^ , since 
then the contribution from all other singularities (for example nucleon self-energy 
cut or anti-nucleon pole) are suppressed at least as (A^/Md) 2 ^. Since the mean 
value of the energy shift is small (A^/AL rm r) oc O(10~ 2 )), Eq.f2|) provide deuteron 
Drell-Yan cross-section with accuracy up to terms of order (A^/Md) 2 cx O(10 -4 ). 



3. Numerical results 

Using the Gluck-Reya-Voght (GRV) parameterization^ for the c pN and solution of 
the Bethe-Salpeter equation with the Graz-II separable kernel^ as / N / D (AfD, p) we 
make numerical calculation fo the ratio <7 pD /cr pp . It is important to note that the full 
GRV parameterization is obtained within the approximation er pD = (cr pp + cr pn )/2, 
therefore it already implies isospin asymmetry of the sea-quark distribution in the 
nucleon in order to describe the DY data. To check influence of the relativistic nu- 
clear effects on the deviations of the ratio cr pD /<r pp from unity we use the GRV 
parameterization without quark-sea asymmetry (u/d = 1). The results of the cal- 
culation are presented in Fig[TJ The dashed curve represents the calculation within 
the approximation with cr pD = (cr pp + er pn )/2, the solid curve represent calculation 
with the four-dimensional Fermi- motion taken into account. 

It is clear from the comparison with the new FNAL data^, that the relativistic 
kinematics can reproduce the deviation of the cross-section ration without the flavor 
asymmetry of the nucleon sea. It is important to note, that the A-dependence of the 
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deviation is denned by the coefficient A^/Ma in front of the cross-section derivative 
same as in case of DIS. 

Thus, we can conclude that observed flavor asymmetry results from the rela- 
tivistic kinematics of bound nucleons; the bound nucleon structure change have the 
same nature in the DIS and DY processes; the effect has the same ^-dependence 
as EMC-effect. A more detailed study of the ^-dependence of the effect and high- a; 
data can shed more light on the nature of the sea-quark flavour asymmetry in the 
nucleon. 
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Fig. 1. Calculations of the DY cross-section ratio for the deuteron and proton compared with 
FNAL E866 (2001) data. The dashed curve is the calculation within the approximation cr pD = 
(<r pp + <r pn )/2, the solid curve is the calculation with the four-dimensional Fermi-motion taken into 
account. Both of the calculation use the GRV parameterization with sea-quark flavour symmetry 
preserved (u/d = 1). 
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